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Dynamics of the neck degree of freedom during fusioning process between heavy ions 
is studied. Time scales of the three degrees of freedom (the relative distance, the neck 
and the mass-asymmetry) are studied, showing an early equilibration of the neck. This 
means that a di-nucleus formed by the incident combination of ions quickly forms a 
mono-nucleus with a superdeformation during the fusion process and that the other two 
degrees of freedom have to be solved in a coupled way. A brief introduction of Langevin 
approach and dissipation-fluctuation dynamics is also given and of the application to 
the synthesis of the superheavy elements. 



1. Motivation : Historical Background 

It is well known that the collective motions in excited nuclei are governed by the 
potential landscape and dissipation. When the system has to cross a barrier, the 
fluctuation, associated with the dissipation according to the Dissipation-Fluctuation 
theorem, play an essential role. 

For multi-dimensional problems, the Langevin equation appears to be easier to 
solve numerically than the equivalent Fokker-Planck equation also used for the anal- 
ysis of heavy- ion collisions, say, of the fast-fission process^. The Langevin equation 
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was first used in nuclear piiysics for a dynamical description of the fission processPl 
Combined with particle emission, it was applied to analyse the anomalous mul- 
tiplicities of pre-scission neutrons,^ which, together with the total kinetic energy 
of fission fragments, supports a strong friction of the one-body type, (Wall-and- 
Windo^^) rather than the two-body viscosity.^ 

More recently, this Langevin formalism has been used to study the fusion of 
heavy nuclei in order to propose an explanation of the fusion hindranc^ which has 
been experimentally known to exist in fusion of massive systems without theoret- 
ical explanation. The DNC (Di-Nucleus Configuration) formed by the contact of 
two ions of the incident channels has an extremely large deformation and then is 
located outside of the conditional saddle point in LDM (Liquid Drop Model), as 
is shown in Fig. 1. The system has then to cross two barriers to fuse. In order to 
calculate probability of the hindered fusion, we use a Two-Step Model and apply 
it to the synthesis of the superheavy element s.^^^ 




Fig. 1. A schematic illustration of LDM energy surface for heavy systems. The x-axis is the 
distance between two centers and the y-axis is the mass-asymmetry. Two examples of incident 
combinations of ions are indicated for the hot and cold fusion paths, respectively. 



The over-passing probability of the inn er barrier is dominated by the diffusion 
of the distance parameter between two ions !-*^^l-^-^l-^^l-*^^l This is briefly recapitulated 
in section 2 with a simplified model. In section 3, time development of the neck 
motion is analysed with the Smolchowski equation during the fusion process. 

2. Brief Reminder of One-Dimensional Parabolic Barrier 

Reducing the inner barrier to an inverted parabola and assuming that the transport 
coefficients are constant near the barrier, the problem is amenable to a simple 
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analytic explanation of the origin of the fusion hindrance. Here, we will recapitulate 
one-dimensional case and refer to Rem for the N-dimensional case with a coupled 
Langevin equation. 

The equation to be solved is given as follows. 
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where fi and uo denote the inertia mass and the frequency of the parabola, respec- 
tively. The reduced friction P = j/fi is defined with the friction coefficient 7. R 
denotes a Langevin force associated to the friction 7, which is assumed to sat- 
isfy the dissipation-ffuctuation theorem, and to be a Markovian with a Gaussian 
distribution. The probability of passing over the barrier, which we call formation 
probability of the compound nucleus, is simply given by an error function, whose 
argument is expressed by the average trajectory < q{t) > and its variance, 

dq 1 _f {q-<q(t)>f 
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For a time long enough, the probability converges to a finite value. 
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where B = /ic<;^(7o/2, the saddle point height measured from the initial point qo, 
while K = Pq/{2/Jj). x denotes the critical parameter l3/(2uo). The probability be- 
comes when the initial kinetic energy K = {x ^ \J x^ + 1)^^, which we call an 
effective barrier 5e// for the case of dissipative dynamics. This simply explains the 
origin of the hindrance. As we discussed elsewhere,"'"^ ^ the distribution of po is 
expected to be a Boltzmann distribution, and then, an averaging over the initial 
momentum po that is thermally distributed with an average value equal to zero, 
gives an extremely simple expression for the formation probability. 
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As is clearly seen in Eq. ([5]), even if we give a larger incident kinetic energy, the 
formation and then the fusion probability increase very slowly through the increase 
of the temperature of the system, which appears in agreement with the experiment .1^ 
Before proceeding to a discussion on the motion of the neck degree of freedom, 
it is meaningful to have a close look at the time evolution of the fusion. We analyse 
time-developments of the trajectory, the formation probability and the current over 
the saddle by the use of the analytic solution. The results are shown in the first, 
the second, and the bottom rows of Fig. 2, respectively. 
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Fig. 2. Average trajectory, over-passing probability and current at the top of the barrier as a 
function of time for four regimes, K = (first column), K = Be///2 (second column), K = B^ff 
(third column) and K = ^B^ff (last column). For each case, two temperatures were chosen, 
T = B/b (solid line) and T = B/2 (dashed line). Note that each column has a different time scale. 



Our case corresponds to the leftmost column, for incident kinetic energy = 0. 
Firstly, the average trajectory never passes over the saddle, but retreats back. Sec- 
ondly, the formation probability becomes saturated as time goes. Finally, the cur- 
rent has a peak structure around several fi/MeV, which shows that the probability 
current at the saddle starts slowly and then terminates gradually. This means that 
the formation of the compound nucleus is not due to a dynamical motion, but due 
to the fluctuation, i.e., due to a tail part of the Gaussian distribution around the 
average trajectory. The time scale of the radial fusion is important for the discussion 
of the neck degree of freedom in the next section. 



3. From Di-Nucleus to Mono-Nucleus : Filling-in of the Neck Cleft 

For a description of fusion processes between heavy ions, there are at least three 
parameters or variables. In two-center parameterization,-'-^ they are distance be- 
tween two mass centers, mass-asymmetry, and the neck correction. Since the neck 
degree of freedom is weakly coupled to the others, it is meaningful to analyse its 
time evolution separately. The LDM potential for the symmetric incident systems 
turns out to be approximately linear in the neck parameter. To analyse the time 
evolution of the neck parameter, starting at e = 1.0 or around, a Langevin equation 
is solved. It appears that the average value of the neck parameter changes very 
quickly, far quicker than the radial fusion for most systems including very heavy 
ones.-*-^ This is due to the action of the linear driving force in the neck e, while the 
radial fusion is governed by diffusion. Thus, it is inferred that the neck degree of 
freedom is in the thermal equilibrium during the fusion. 

Next, in order to know how the distribution reaches the equihbrium, we try 
to obtain a time-dependent distribution function of the neck, starting from the 
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delta- function at eo = 1.0, i.e., at the initial DNC. The Smoluchowski equation is 
solved, since we know that the momentum space can be approximated to be in an 
equilibrium, due to a very small inertia mass."*"^ Then, with a linear potential, the 
equation to be solved is as follows. 



dN 
~dt 



D- 
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where the diffusion coefficient is D = T and the drift one C = f /'^^ f being the 
slope parameter calculated with LDMI^ ^(^) = The friction coefficient 7 is 
calculated with the usual one body model. For simplicity, we take the range of the 
variable e to be [0.0, 00], instead of the realistic [0.0, 1.0] (in this case, a little more 
complicated expression has been obtained, but the results are essentially the same 
as the present case.*^). The boundary condition at e = 0.0 is reflective. With the 
initial and the above boundary conditions, the solution is obtained as follows 
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For long times, this expression becomes a Boltzmann distribution. In Fig. 3, the 
time dependence is shown by distributions at various times after the contact, for 
the case of ^°°Mo+^°°Mo system. 



Time evolution 

of the neck distribution 




Fig. 3. Time evolution of the neck distribution function is shown with examples at various times 
for Mo+-'^ Mo system, for which typical values of the parameters are D = T/j=l/8 and 
C = //7=20/8 in the unit of MeV/h. The time unit is h/MeV. 



Apparently, the Boltzmann distribution in the coordinate space is established 
at several tenths of h/MeV. It is worth noticing that this time scale is far shorter 
than that of the 1-dimensional radial fusion discussed in the previous section. This 
means that before the radial motion for fusion starts, the neck cleft is ffiled in, i.e.. 
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the initial DNC becomes a superdeformed mono- nucleus. During fusioning motion, 
we can approximately take e to be 0.0, because it is the most probable value of the 
Boltzmann distribution obtained above. 

A similar analysis has been made for the mass-asymmetry,-*-^ which turned out 
that the time scale is the same order with that of the radial fusion, and thus, the 
two degrees of freedom have to be solved in a coupled wayP 
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